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1 Introduction 



Ever since its invention super symmetry has been an important subject in high-energy 
physics beyond the standard model. It is considered to be a necessary ingredient to 
bridge the gap between the scale of electroweak symmetry breaking and the much larger 
unification scale. Nowadays, super symmetric theories cover the whole range from su- 
persymmetric classical mechanics jlj , quantum mechanics (21 |S1 , scalar and gauge field 
theories [1] to string- and M-theory [3] . They allow for the construction of low-energy 
effective actions, as for the M = 2 Seiberg-Witten model [£] or the formulation of certain 
duality relations, like in the original Maldacena conjecture for gauge theories with M = A 
extended supersymmetry 0. 

The non-perturbative effects in supersymmetric theories, and in particular, the dynam- 
ical breaking of supersymmetry are a subject of intensive studies. At present time the 
lattice formulation is the only tool for systematic investigations of such effects, and lat- 
tice simulations provide the means of doing reliable calculations in the strong-coupling 
regime or near a phase transition point. After the pioneering work of Dondi and Nico- 
lai [S] there has been an ongoing effort into formulating, understanding and simulating 
supersymmetric theories on the lattice [HlES^DEl- Recent lattice results, e.g. on the 
breaking of supersymmetry, have been obtained in |131 1141 [T^ . 

A commonly accepted guiding principle in any good lattice calculation is to build in as 
many of the symmetries of the continuum model as possible, such that the lattice results 
respect these symmetries identically. However, often these are conflicting requirements 
and not all symmetries can be incorporated on the lattice. This in turn introduces subtle 
lattice artifacts into the formalism, which one may not get rid of in the continuum limit. 
For example, lattice regularizations of supersymmetric theories generically break large 
parts of supersymmetry, and it is a nontrivial problem to recover supersymmetry in 
the continuum limit. However, there are discretizations with highly nonlocal derivative 
operators, for which supersymmetry is manifestly realized [HI Alternatively, for two- 
dimensional models one can discretize only space (time remains continuous) such that a 
subalgebra of the N = 1 supersymmetry algebra, 

{Qa,Q/3}=2(7'^7%/3^M 

remains intact [HIEIE]- That subalgebra then determines the spectral properties of 
the super-Hamiltonian H. The fermion doubling for naive lattice derivatives jl7| 
is another apparently unrelated notorious example of such lattice artifacts. For bosons 
there is no such problem. However, if we try to preserve part of supersymmetry then 
the fermionic mirror states lead to doublers in the bosonic sector as well. 

In this paper we study continuum and lattice versions of two-dimensional Wess-Zumino 
(WZ) models. Similar to the original four-dimensional theory pjij, these models contain 
scalar and fermion fields coupled by a Yukawa term. A particular version possesses 
M = 2 supersymmetry and has been the subject of analytic HOI |^ and numerical [22] 
studies. 
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In section 2 we consider the off-shell formulation for a general class of continuum mod- 
els and derive the supersymmetry transformations and Noether currents. Particular 
emphasis is put on the form of the central charges j23j . 

In section 3 we turn to the lattice version of the models. We show that for real and 
antisymmetric lattice derivatives the M = 1 algebra can be represented on free fields. 
The local left- and right-derivatives are not antisymmetric and the anticommutator of 
the corresponding supercharges does not yield the discretized Hamiltonian for the free 
model. If we insist that supersymmetry is realised on free fields without fermion and 
boson doubling then we must allow for nonlocal derivatives on the lattice. One particular 
such derivative, the SLAG operator, is introduced in this section. The numerical results 
for this operator concerning supersymmetry in lower-dimensional systems are in excellent 
agreement with continuum results. In section 4 we show how to derive the models with 
M = 1 and J\f = 2 supersymmetry on a spatial lattice by a suitable dimensional reduction 
of a high-dimensional Euclidean Dirac operator. In the process of reduction the Dirac 
matrices and coordinates turn into Majorana spinors and scalar fields on the lattice. 
We count and construct the normalisable eigenstates of H with zero energy both in the 
weak and strong-coupling limits. In particular we find that the M = 2 models with </>^'^ 
interaction admit such states if N is the number of spatial lattice sites. 

In section 4 we bridge the gap between strong- and weak-coupling regimes for models 
with M = 1 and N = 2 supersymmetry with the help of powerful methods from operator 
theory. Using a theorem by Kato we prove that the zero modes in the strong-coupling 
limit survive for intermediate couplings as long as the coupling constant of the leading 
term in the potential does not vanish. We comment on what we expect to happen in the 
continuum limit of the M = 2 models, where only q of the q^ zero modes survive [21]. We 
also comment on recent lattice simulations of the two dimensional Wess-Zumino model 
by Beccaria et al. |25j . Some technical details concerning the nonlocal SLAG operator 
and the proof that the transition from strong to intermediate couplings is governed by 
a relative compact perturbation are relegated to the appendix. 

2 Wess-Zumino Models in 1 + 1 Dimensions 

In the off-shell formulation two-dimensional parity invariant Wess-Zumino models con- 
tain a set of, say d, triples, each containing a real scalar (j), Majorana spinor -0 and 
auxiliary field F. In a Majorana representation for the Glifford algebra 

{7^7'^} = 27?'^^ with /^Mt/ = ^A«, ,y = diag(l,-l), (1) 
the Majorana spinors are real. 

The supersymmetry algebra is spanned by M Hermitian spinorial supercharges Q'^^\ 
I = 1, . . . , AA, by the Hermitian two-momentum and by the (anti-)symmetric matrix 
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of Hermitian central charges Zg"^ {^a) ^^'^ form 

{Q'a\Qf} = 2 {S'-^Kp + i5a/3^A^ + i^a/^^-") > 7* = , (2) 

with spinor index a = 1,2. 

In component fields the Lagrangian of the models with M = 1 supersymmetry reads 

where the superpotential W depends on the dimensionless scalar fields (j)^,... ,4>'^. We 
denoted the derivative of W with respect to by W,a and employed the Einstein sum- 
mation convention. For Wess-Zumino models the target spaces are R"^ with Euclidean 
metric Sab- 

Now we consider the most general linear off-shell supersymmetry transformation of the 
fields. Since {(/)"', i/j"' , F'^) have mass dimensions (0,^,1) respectively, such transforma- 
tions have the form |26j 

6,i;'' = \${B(t)Te + {CFYe, (4) 
<5,F» = i-e${D,l,f, 

where, for example, {AiJjY = A^-^il)^. The constant matrices A,B,C,D must be real for 
the supersymmetry variations to be Hermitian fields. The requirement that C transforms 
into a divergence implies the following algebraic relations for these matrices and the real 
symmetric matrix W" = {W,ab ), 

A + B^ = 0, D + C^ = 0, (5) 

A^W" + W"C = 0, W'A'^ + CW" = 0. (6) 

It follows that 

5C = ed^V + A with A = -^W,abc {eA'^^i^'') (^V') • 

Free models have quadratic superpotentials and A is identically zero. For interacting 
models we may exploit the Fierz identity 

to prove that A vanishes, provided 

W"A = A^W" (7) 

holds true. Then the action is left invariant by the transformations (0]) and the corre- 
sponding conserved Noether current reads 

J^' = {^^'(t)--i^e^''^^(t)T{A'4:Y --liCW'TY^". W = {dW / dtt)"") . (8) 
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In what follows, employing ©, we express the matrices B and D in terms of A and C. 
We consider M supersymmetries @ with matrices {Aj, Cj) and denote the corresponding 
supersymmetry transformations by 5i^\ For all pairs {Ai,Cj) the conditions © and 
(O must hold for the Lagrangian to be invariant. These conditions severely restrict the 
form of the superpotential W. We also demand that two supersymmetry transformations 
close on translations (later we shall comment on the possibility of central charges) 

[5(f),<5(f]cl> = 2i<5"(6-27^6i)a,cl>, (9) 

and this puts further restrictions on the matrices. For the scalar and the auxiliary field 
the condition @ read 

AjA^ + AjAJ = CfCj + CjCi = 26ijl and AjCj - AjCj = 0. (10) 

In particular all matrices are orthogonal, such that the two conditions in © coincide. 
Actually, the last relation implies that the algebra © is realized on the Majorana fields 
as well. 

The transformation di^^ is generated by the Noether charge corresponding to Jj in ©, 
Q(^) = Jdx ((vr - cP'-f,)aiAJ^pr - i(C,W^')a7V) , vr'^ = (H) 
where we have set {d(j)'^/dx) = (j)' . 

Canonical structure: The canonical structure is more transparent in the on-shell for- 
mulation. This is obtained from the off-shell one by replacing Fa by W,a- The nontrivial 
equal time (anti)commutators between the scalar fields, their conjugated momentum 
fields tt"" = (f)"" and the Majorana fields read 

{C(x), V'^(y)} = 5^p5''H{x - y) and [r{x),7:\y)] = id^'dix - y). (12) 

The Hamiltonian is the Legendre transform of the Lagrangian, 

H = j dxH, n = \Ti-^ + \(i)' -4)' + \w' -w' + (13) 

where, for example, vr-vr = vTaVr". We have introduced the Hermitian Dirac-Hamiltonian 

{h^)ab = -h*dJab + 7°W^,afe = {K)ab + 7°W^,a6 ■ (14) 

The action is invariant under spacetime translations generated by Noether charges 

Po = n and Pi = J dx (vr • (/>' + i^7V) , (15) 

and under supersymmetry transformations (0]) generated by the above supercharges Q^^\ 
By using the relations H6I1U() one proves that the Q^^^ satisfy the super-algebra @ with 
central charges 

Z'/ = and Z's^ = -j^'- {AjCj)W', (16) 
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where we have neglected ambiguous surface terms containing the Majorana fields only. 
Note, that the integrand is a total derivative, since the integrability conditions for the 
existence of a potential U{(j){x)) with 

</.' • {AjCj)W' = ^ = [/' . 

is that AjCjW" is a symmetric matrix. But this follows from the condition 
In most explicit calculations we choose the Majorana representation 

7° = <72, 7"^ = io"3 and 7* = ■j^j'^ = -a\ (17) 

such that the superalgebra takes the simple form 

{q(^),q('^)} = 2{Bb'-^ -Z's'). (18) 
{qS^\q(^)} = 2 (Pi^^-^ + Zi^) . 

M = \ supersymmetry: There is always at least one solution to the constraints (|5|(il7|) 
and ()1U|) for an arbitrary superpotential VF, namely 

Ax = -Bi = -Ci = Di = \. (19) 

Solving for the auxiliary field. Fa = W,a, the on-shell transformations take the form 

6P(j) = eip, dPip = - W') e, (20) 

and the corresponding supercharge reads 

= J dx {tt- (j)'-f^ + iVF'7°) • iP . (21) 

For vanishing spinors the only non-trivial central charge is 

d.i5:. (22) 

dx 

J\f = 2 extended supersymmetry: We assume that the model © admits a second 
supersymmetry besides the solution (fT^ . The conditions (fTU]) imply 

A2 = -C2 = I, I = -I^, I^ = -l. (23) 

The matrix / defines a complex structure and exists for all target spaces IR*^ with even 
dimension d. The conditions in ® and ((7j) on the superpotential both reduce to 

IW" + W"I = 0, (24) 
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which means that the superpotential is a harmonic function of the scalar fields, in agree- 
ment with the general analysis in On-shell, the second supersymmetry has the 
form 

<5f V = e^V', <5f V = {i^I(t> - IW) e, (25) 
and is generated by the Noether-supercharge 

q(2) = j dx {tt- 0'7, - iVF'7°) • (ItP). (26) 

For vanishing spinor fields the central charges read 

Z'/ = and {Z's') = ^3 y dx ^ - C71 y dx ^, (27) 

where U is the imaginary part of the analytic function ^(0-*^ +i(/)2) = W + iU with real 
part W. 

For the models with N = 2 supersymmetry there exists a concise formulation in which 
two real scalar s are combined to a complex scalar, and two Major ana spinor s are com- 
bined to a Dirac spinor. For example, for the target space we set 

(/.= i=((A^+i,/.2), V' = ^(V''+i7*V''). (28) 

The harmonic superpotential is the real part of a holomorphic function, 

Ty(0,0) = F(</))+F(0), (29) 

and the on-shell Lagrangian takes the form 

£ = d^c/jdf'cl)^ + iip0^ - - F"i)P+i} - F"i)P.i), (30) 

where F' is the derivative of F with respect to the complex field ^ and we have introduced 
the chiral projectors 

P± = i(l+7,). (31) 

Along with the real scalar fields one combines the corresponding conjugate momentum 
fields to a complex momentum, tt = (vr^— i7r2)/\/2, such that 

[(/)(x),7r(y)] = i5(x - y) and {V'a, V^/j} = ^^a/?- (32) 
The complex supercharge takes the form 

Q = \ (q(^) + i7*Q^^^) = (tt - + iFV) P+ij + (vf + 0' + iF'7°) P.^. (33) 
and satisfies the anticommutation relations 

{Q,Q} = and {Q,Q] = f + -f.Zf - Zf . (34) 
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Higher supersymmetries: Next we show that with the absence of central charges 
there is no third hnear off-shell super symmetry besides H2U|) and ()25() . To be compat- 
ible with the first transformation in (|2U() . the orthogonal matrices ^3 and C3 must be 
antisymmetric and of opposite sign. The conditions H1U|) between the second and third 
supersymmetry imply 



which is impossible for orthogonal matrices I and A^. We conclude that the models (jSJ 
admit at most two linear off-shell supersymmetries. 

Let us mention that, if we allow for central charges in the superalgebra, there exist 
further supersymmetries. But the corresponding models are massive free models. They 
can be derived by a dimensional reduction of the free M = 2 model in 4 dimensions. 

3 Lattice Formulations of Wess-Zumino Models 

As ultraviolet-cutoff we discretize space, introduce a spatial lattice with equidistant 
sites and choose periodic boundary conditions. The time is kept continuous such that 
time translations remain symmetries generated by the Hamiltonian. Following [2] we try 
to preserve at least that subalgebra of © which involves H. 

The fields of the supersymmetric model in the Hamiltonian formulation are discretized 
as follows. 



where the lattice spacing has been set to one. On a space-lattice the derivative becomes 
a difference operator the particular choice of which is left open for the moment being. 
We define the lattice Hamiltonian as square of the discretized supercharge Qi. For 
interacting theories it consists of the discretized Hamiltonian of the continuum theory 
plus a lattice counterpart of the central charge. 

On-shell the N = 1 model contains d G {1,2, . . . } Hermitian scalar fields (l)"'{n) and d 
Major ana spinors ifj'^in) on N lattice sites (n = 1, . . . , N). The fields obey the non-trivial 
canonical (anti-)commutation relations 



[4>''{n),TT\n')]=\5''H{n,n') and {C(n), '0^(n')} = 5"X/35(n, n'). (36) 



We choose a Majorana representation such that the are Hermitian two component 
spinors. 

When we put the supercharge on a space-lattice, we must choose the lattice derivative 
in the term 



[7,^3] ={1,^31 = 0, 



{4>''{x),7^%x),r{x)) 



{r{n),7^%n),r{n)) 



n = l,...,iV, 



(35) 
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in (jllj) . Since we do not want to specify d at this point we make the general ansatz for 
the Hermitian Dirac-Hamiltonian 

hl = -i5abl^^Q^ , With dd^ = d^d = -A (37) 

and a real d with correct continuum limit, d must be real, since it should map Majorana 
spinors into Majorana spinors. Let us define its symmetric and antisymmetric parts 

9s = i(9 + at), d^ = \{d-d^) with [5a,9s] = 0, dl-dl = ^. (38) 

The last two properties follow from our assumption [9,5^] = in ()37() . Since 

/,o*P-i^,aA-7°<9s, (39) 

chirality is preserved for massless fermions if d = is antisymmetric, in which case 
hp = — i7*5A- Thus, if d is antisymmetric and local then, according to some long- 
standing no-go theorems there is fermion doubling. There are many such theorems, and 
we mention only two, one due to Nielsen and Ninomiya |T7] and a later elaboration due 
to Friedan ^Hl- No-go theorems are notorious in that people find a way around them, 
and following Friedans work, Liischer ^| and others did so. Below we circumvent the 
no-go theorems by using a nonlocal and antisymmetric derivative. 

However, most lattice derivative are not antisymmetric in which case hp contains a 
momentum dependent mass term —j^dg. Such a chirality violating term has been in- 
troduced by Wilson to raise the masses of the unwanted doublers to values of order 
of the cutoff, thereby decoupling them from continuum physics. 

As discretized supercharge (|^T|) we take 

qW = (vr, + i(</), h^i;) + [{W , 7V). (40) 
A careful calculation yields the following anticommutation relations, 

UQa\Q^^^} = {flXp - i{.lXp{W\d^<t>) - 5^p{W\dsct>) (41) 
with energy and momentum 

2Po = {7T,n)-{(t>.M) + {W\W') + {il,,hpij), 

2Pi = 2(aA0,7r)-(V,7*/i?V'), /ip = /i°+7°^^"- (42) 

To arrive at these results one uses the identity 

(7r,a0) +i(at^i,Vi) = (a</.,^) -i(Vi,5tV'i), 

which holds for any real difference operator d. The superalgebra can be rewritten as 

= f + \^,{W\ 5a0) - l\W', 9s0). (43) 
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The last term is absent in the superalgebra © and breaks Lorentz covariance explicitly. 
This lattice artifact originates in the Wilson term —^^ds in H39() . This term must vanish 
in the continuum limit. One may wonder whether there exist other improvement terms 
we could add to a local — i7*9A in order to avoid the fermion doubling. However, since 
for Major ana fermions the terms 

(V',9sV'), i^,'y^ds'ip), (V',7*5s^) 

are constant or zero, all terms but j^ds do not show up in the right hand side of ()42() 
and we obtain the same result as if we had chosen hp = —i'y^:dA- Hence only the Wilson 
term ~ 7^ 9s can be used to avoid the fermion doubling. This argument does not apply 
to theories with several Majorana fermions and in particular to models with extended 
supersymmetry. 

Models with J\f = 2 supersymmetry contain the second supercharge in (|26|) . the lattice 
version of which reads 

= (tt, I^) + hp^) - \{W\ ^^I^), (44) 

and satisfies the same anticommutation relations as Q*-^^ up to a sign change of the last 
two terms in ()48|) . The anticommutator of two lattice charges reads 

1 {Q(^), Q(-^)} = 5''f + ij^Z's' + Zi\ (45) 

where the 'would-be' central charges 

Z's' = ^i'iW, 9^0) - (ai)"(f/', dA<P) (46) 

approach the central charges H27() of the continuum model. To arrive at H45() one needs 
the harmonicity of the superpotential which in turn implies the existence of a function 
f7((/>) with IW = U' , and this function enters the central charges. However, since the 
Leibniz rule never holds on the lattice, the integrands W ■ dx(t> and U' ■ d^fp in (|46|) are 
not just total derivatives as in the continuum and as a consequence the terms Z^"^ are 
not central to the algebra. The annoying terms 

= -(a3)"7°(^',9s</') + (^i)'-' {7'{U',ds<P) - i(vr,/5s0) - |(V,/5sV)) (47) 
in 1)451) are pure lattice artifacts and vanish for antisymmetric lattice derivatives. 

Free Wess-Zumino model (AA = 1): For simplicity we consider the free model with 
scalars of equal mass. The superpotential reads W = ■^ttkI) dCf)^ and. with. W^' — ttkI) the 
'would-be' central charge vanishes, 

(Ty',9A</.) =m(</.,5A</.) = 0. (48) 

As Hamiltonian we choose the square of the supercharges, 

H = \{Qi,Qi] = \{Q2,Q2] = P^-m{(l),ds4>), 

2Po = (vr,^) + (,/.,(- A + m2)0) + (V',/ipV'), (49) 
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where the Dirac-Hamiltonian for the non-interacting model is just 

hp = -ij^dA + 7°(m - ds) with hj = (-A + - 2771^3)12 (50) 

and — A = dd^ . For antisymmetric derivatives the pure lattice artifacts containing dg 
vanish and with 2Pi = {Qi, Q2} we obtain the familiar algebra 

{Qa,Qp} = 2{j^jXf3P^, [Qa,Pi,]=0, [Po,Pl]=0. (51) 

We conclude that the = 1 superalgebra in 1 + 1 dimensions can be represented as a 
free Wess-Zumino model on a space lattice. 

3.1 Lattice Derivatives 

At this point some words about lattice derivatives are in order. At first instance one 
may think that the local right- and left derivatives 

(aR,/)(n) =/(n + l)-/(n) and (dMn) = f (n) - f (n - 1) (52) 

are ideal candidates for a lattice derivative. With respect to the ^2-scalar product of two 
lattice functions, 

N 

{f,9) = Y.f{n)g{n), (53) 

n=l 

the adjoint of the left-derivative is minus the right-derivative, dl = — (9r. Both derivatives 
share the property that (1, 9r/) = (1, dj^f) = 0. But the corresponding momenta Pl = 
— i^L and = — ISr are not Hermitian and possess complex eigenvalues, 

XkiPn) = hiP0 = 2(^P>'/^sm^, with pk = 2Trk/N, and k = l,...,N. 

If we insist on a Hermitian momentum we could choose the antisymmetric derivative 
operator 

5r+i. = ^(5r + 50 = -Cl (54) 
which is used in many lattice calculations. The N real eigenvalues of Pr+l read 

Afe(pR+L) = smpk = Re {XkipR)) , 

and waves with the shortest wavelength, that is with at the boundary of the first 
Brillouin zone, are zero modes of 9r+l- Hence, by trying to preserve the hermiticity of 
p in this naive way immediately introduces spurious zero modes that are responsible for 
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the fermion doubling problem. 



• , • • o 
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Aa:(j5r+l), ReXkipB 
Im Xkipn) 
Aa:(j5slac) 



A third alternative for the lattice mo- 
mentum is the Hermitian and nonlocal 
SLAC operator Pslac = — i<9sLAc intro- 
duced by Drell, Weinstein and Yankielow- 
icz with real eigenvalues pk- This op- 
erator has no spurious mirror states in the 
first Brillouin zone. The eigenvalues of 
PsLAc are eigenvalues of the momentum 
operator in the continuum and the differ- 
ence between the lattice and continuum 
results are minimized. In the figure on the 
I left we have plotted the eigenvalues of the 
various lattice operators. The real parts 
of the eigenvalues of and Pl are just 
the eigenvalues of Pk+l- The eigenvalues 
of Pn+L are twofold degenerate. The SLAC 
operator has the same dispersion relation 
as the momentum in the continuum. 



Besides d^, c^l, ^r+l and c?slac there are many other local and nonlocal candidates for lat- 
tice derivatives with the correct naive continuum limit. However, it is easy to see that no 
linear difference operator will obey the Leibniz rule. Many problems in supersymmetric 
lattice theories are exactly due to this fact, see 0. 

In order to better understand the dependency of the spectrum and doubling phenomenon 
on the lattice derivative we consider the following one-parameter interpolating family of 
ultra-local difference operators 



da = ^(1 + a)d^ + i(l - a)(9L = ds + 5a, 
with symmetric and antisymmetric parts 

= ^a{d^ - (9l) = ^ad^di^ and 5a = ^(5r, + 3l) = Sr+l. 



(55) 



(56) 



When the parameter a varies from 1 to —1, then da interpolates between and (9l- 
For a = we obtain the antisymmetric operator 5a in ()52() . 

The 2A^ eigenvalues of the Hermitian Dirac-Hamiltonian (|5fl|) depend on the deformation 
parameter as follows. 



Afc(a) = A7v_fc(a) = m? + Aa{a+m) sixiWpk) + (l-a^) sin^(pfe), (57) 
where pk = 2'Kk/N and k runs from to A^— 1. For the extreme cases a = 0, 1 we obtain 



12 



the eigenvalues 
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Afc(O) = ±^m^ + sm^ (pk) 
with multiphcity 4 and 



Afc(l) = ±^Jm^ + 4(1 + m) sm\^pk) 

with multiphcity 2. This should be com- 
pared with the eigenvalues on the contin- 
uous interval of 'length' N, 



(58) 



with multiplicity 2. One can show that for 
a greater then a+ or less then a_, where 



4a± = ± ( Vm^ + S =F m] , 



all eigenvalues have the same multiplicity as in the continuum. In particular, for massless 
fermions there are no doublers for > 1/2. However, for a G a+] some eigenvalues 
have multiplicity four. In the above figure we have plotted the positive eigenvalues of hp 
for a = 0, 1, a+. For comparison we have depicted the positive eigenvalues of hp for the 
nonlocal SLAG derivative 



(a 



SLAc)n=in' 



tt/N 



sin (7r(n — n')/N^ 



and (9s 



0. 



(59) 



Despite being nonlocal the SLAG derivative has many advantages as compared to the 
local operators d^,di^ or c?r+l: it is antisymmetric such that for massless fermions chiral 
symmetry is preserved. By construction the 2A^ real eigenvalues of hp = —i^/^dsi^Ac+l^fn 
are identical to the 2A^ lowest eigenvalues of the continuum operator on the interval of 
'length' N, For this reason ds^Ac has been called ideal lattice operator in the 

literature. We do not expect that unwanted nonlocal counterterms ^30^ are required 
for the two-dimensional supersymmetric Wess-Zumino models. This is certainly the 
case for the finite models with extended supersymmetry. For the model with J\f = 1 
supersymmetry the same should be true since it does not contain gauge fields which 
couple to high momentum modes at the edge of the Brillouin zone. Indeed, in is has 
been claimed that (?slac approaches an ultra-local operator when N tends to infinity, 
except for a border matrix. In the appendix we give a detailed analysis of this interesting 
operator. 
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3.2 On the Quality of Lattice Derivatives in Supersymmetric QM 



It is enlightening to retreat to quantum-mechanical systems and study the supercharges 

Q = (^^t Q ) , with A = d + W, = d^ + W, (60) 

and in particular the quality of lattice approximations for different lattice derivatives d 
in A. The supercharge squares to 

with isospectral discretized Schrodinger operators 

AA^ = dd^ +dW + Wd^ + W'^ 

A^A = d^d + d^W + Wd + W'^. (62) 

They have identical spectra, up to possible zero modes. If the Leibniz rule held on the 
lattice, if d was antisymmetric and if we could replace dW by W' + Wd, then we would 
find the super-Hamiltonian of supersymmetric quantum mechanics in the continuum, 

_ (dd^ + w' + w^ \ 

d^d-w' + w^)- ^^'^^ 

The difference between and H is the analog of the last two terms in H43|l and the 
difference in their spectra is a good measure for the suitability of the chosen lattice 
derivative as regards supersymmetry and the speed with which the continuum limit is 
approached. In the following figure we have plotted the eigenvalues of and H for 
d = dsLAc, denoted by QIlac ^"^^ -f^SLAC and for d = d^,, denoted by Ql^-^ and -ffnaiv 
We took the superpotential W = Xx^ which gives rise to the supersymmetric anharmonic 
oscillator. 

The lowest 57 eigenvalues of and H are almost identical for the SLAC derivatives and 
the lowest 90 eigenvalues of -ffsLAC agree with the exact values (calculated on a much 
finer grid). These results clearly demonstrate the high precision of the SLAC derivative 
in low-dimensional supersymmetric systems. It does not matter whether we discretise 
the supercharge or the super-Hamiltonian as long as we choose the SLAC derivative. 
After this detour to quantum mechanics we now return to supersymmetric field theories. 



4 From the Dirac Operator to the Lattice M = I WZ Model 

In this section we relate the supercharges and Hamiltonians of two-dimensional Wess- 
Zumino models on a spatial lattice to suitable Dirac operators. We shall use the results 
in IS2] on the (extended) supersymmetries of \Ip in arbitrary dimensions, specialized to 
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Figure 1: Eigenvalues of and H for the SLAC derivative and the right-derivative for 
N = 180 lattice points, length L = 30 and A = 1. 



flat space and perform a dimensional reduction such that the supercharges of the lattice 
models can be related to the reduced ilp. 

To see how Dirac operators relate to multi-dimensional supersymmetric matrix-Schrodin- 
ger operators we generalize the reduction of the two-dimensional Dirac equation to the 
Nicolai-Witten operator to higher dimensions. For that purpose we dimensionally reduce 
the Euclidean operator 

i^ = ir^i?^, D^ = d^ + iA^, r^^ = r^, /x = i,...,2iv, (64) 

from the products of cylinders, 

M = R X ■ ■ ■ X R X X ■ ■ ■ X 5\ = X (65) 

A''-times A''-times 

to the factor R,-^. We take coordinates on R'^ and 6^, . . . , 9^ as coordinates 

on the torus , respectively. We dimensionally reduce by assuming that the Abelian 
gauge potential is independent of the angles 0". Then the Dirac operator commutes 
with the (angular)momenta —idgn and we may set dgn = 0. 



4.1 Reduction to Models with M =1 

If we further set Ai{x) = ■ ■ ■ = Aj\f{x) = and define n = N+n, then the square of the 
reduced Dirac operator takes the form 

-Ip^ = 2H = pnPn + AnAn-ir''r'^dnArn, where p„ = lA. (66) 

1 ox"- 
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Note that the reduced operator H contains no first-derivative terms. It can be identified 
with the Hamiltonian of a two-dimensional J\f = 1 WZ-model on a spatial lattice, if we 
set 

pfi ) = V2i^{n). (67) 

It follows with (|17() that {ijj^'^^il)) = — iF^F" holds true. If we further assume that the 
non- vanishing components of have the form 

An = -{dU){n) + W'{(p{n)) with = ds - dj,, (68) 

then we find 

- 10 = (vr, y^i) +^{W', ^2) -(</>, aV'2)^ = Qi'^ (69) 
Ao Ai 

with Q^-^^ in H40|l . We conclude that the Hamiltonian reads 

H = Po + {W', dj^ct)) - {W, dsc/)), (70) 

with Pq from (|42() . Thus we have proved that the super-Hamiltonian of the M = 1 
Wess-Zumino model on a lattice with N sites is just the square of the Dirac operator in 
2N dimensions, dimensionally reduced from x to R^. 



4.2 Ground State of the Free Model 



For the massive non-interacting model we have 2W = mclP'. The corresponding Hamil- 
tonian is the sum of two commuting operators, of the bosonic part 

H^ = \{TT,-K) + \{(t),A^<P), A^ = -^ + mds + m^, (71) 

and the fermionic one 

Fp = i(V',/ipV'), /ip = -i7*9A + 7°("i-9s). (72) 

We assume that the parameters are such that is positive. Near the continuum limit 
this is always the case if the physical mass is positive. The ground state wave function(al) 
of the supersymmetric Hamiltonian factorizes, 

^'o = ^'b^f with i^B^B = -E'b^'b and i^p^'p = Sp^'p. 
We choose the field representation for the scalar field, such that 

^W = -7^:|-T and ^^ = m^[4>). (73) 
1 omyn] 
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The bosonic factor ^'b is Gaussian 

^'b = c • exp {-^{(j), A(j))) and = ^tvA. (74) 

Here A is the positive root of the positive and Hermitian in H71|) . For the family of 
operators in (|55jl the trace of A is just half the sum of the positive eigenvalues in (|57j) . 

To find ^I'p we introduce the (two-component) eigenfunctions of hp with positive 
eigenvalues. Since the Hermitian matrix hp is imaginary the cannot be real and we 
have 

hpVk = XkVk hpVk = -XkVk (Afc > 0). (75) 
The eigenvectors are orthogonal with respect to the Hermitian scalar product, 

{vk,Vk')= ^ Vka{n)vk'a{n) = 5kk' and (-y^, ti^,) = 0. (76) 

n,a=l,2 

Now we expand the Majorana spinors in terms of this orthonormal basis, 

N 

= ^{xkVk{n) +xivk{n)^ , where Xk = {vk,'^), xi = {vk,ip) (77) 

k=l 

are one-component complex objects with anticommutation relations 

{Xfc,Xfc'} = and {Xk,x\i] = hk' ■ (78) 
Inserting the expansion (|77j) into Hp yields 

H. = hYl >^k[xixk-Xkxi)- (79) 

fc:Afc>0 

It follows that the ground state of Hp is the Fock vacuum which is annihilated by all 
annihilation operators xi , • • • , Xa^ and has energy 

Ep = -\ (80) 

fc:Afe>0 

Since /i| = I2 A'^ we conclude, that the positive eigenvalues of hp are identical to the 
eigenvalues of A such that 

E = Eji + Ep = 0. 

Since ^'o is normalizable for A > Owe see that the Hamiltonian admits a supersymmetric 
ground state for all choices of the lattice derivative d, provided A is positive. 
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4.3 Ground State for Strong Coupling 



In the extreme case of very strong self-coupling of the scalar field we may neglect the 
derivative term in the supercharge (|69)1 jH] • Then Q and H are the sum of N identical and 
commuting quantum mechanical operators, each defined on a given lattice site. Hence, 
the ground state is a product state, ^'o(</') = ®n^o('/'n)- The operators on a fixed lattice 
site read 

= j^iQ^ + i'2W'{<l,) and H = -^ + W'^ -\xl,ii^2W". (81) 

A normalizable zero-energy state is annihilated by Q^i\ 

M^)=e-''^'^''^^'^^^o, (82) 

where ujq is a constant two-component spinor. It is well-known jSHl and follows at once 
from ()82() that supersymmetry is unbroken if p in 

W = ^X(j)P + 0{(l)P-^), A^O, (83) 

is even and it is broken if p is odd. Note that —iipi'ijj2 is Hermitian and has eigenvalues 
±1 and that for even p the state V'o is normalisable if 

i sign {X)7piTp2u;o = ujq or {ipi - i sign (A)V'2) t^o = 0. 

To summarize, for even p the M = 1 Wess-Zumino model on the spatial lattice has 
always exactly one normalizable zero mode in the strong-coupling limit. For A > this 
product state has the form 



^o(</') = exp (^-5^^VF(0„)^ Oo, 



(V'i(n) -iV'2(n))Oo = 0, Vn. (84) 



In particular, for 0^-models supersymmetry is broken in the strong-coupling limit where- 
as it is unbroken for (/(^-models. 



5 From the Dirac Operator to the Lattice N = 2 WZ Model 

It is known, that on flat spacetime the Euclidean Dirac operator admits two supersym- 
metries if the field strength commutes with an antisymmetric and orthogonal matrix /, 
which defines a complex structure j32j . The two real supercharges 



-^ir^Z)^ and Q2 = ^ 



Qi = — ir^Z)^ and Q2 = — iI^XDu (85) 
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form the superalgebra 

UQuQj} = ^ijH. (86) 
They can be combined to a nilpotent complex charge 

Q = (Qi + iQa) (87) 

and its adjoint , in terms of which the supersymmetry algebra takes the form 

H = Q^}, g2 = Qt2 = and [Q, H] = 0. (88) 
To obtain J\f = 2 lattice models on N sites we consider the Dirac operator on 

M = lRx...xIRxS'^x...xS'^ = R^^ x T^^ (89) 

^ V ■' ^ V ■' 

2Af-times 2Af-times 

in contrast to the 2A'"-dimensional space in (|65|) . Since the field strength commutes with 
the complex structure / it is very convenient to introduce the corresponding complex 
coordinates on M, 

= + ix" = x" + 16''', n = 2N + n, n, n G {1, . . . , 2iV}, (90) 
and fermionic annihilation and creation operators, 

^ 1 ^pn ^ -pn^ ^ V'^" = ^ (r" - IT") with {V'", V'"^"'} = (J™"- (91) 

The condition that F^y commutes with the complex structure / implies the existence of 
a real superpotential xiz,z), such that [221 



/2N ^ \ 



5.1 Reduction to Models with M =2 

Again we perform a dimensional reduction by assuming that x does not depend on the 
compact variables 0^, 

X = x{x\...,x''^) (92) 

and that the angular momenta dgn vanish. In the sector with vanishing angular momenta 
the complex charge simplify to 

Q = e-^Qoe^ where ^0 = 1^^^"^, (93) 

n=l 
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since the complex z"-derivative becomes the real x"-derivative in this sector. This 
dimensional reduced supercharge and its adjoint generate the superalgebra H88|) with 
super symmetric matrix-Schrodinger operator 

H = i{Q, Qt} = - 1 A + i(Vx, Vx) + iAx - V r^X,nm^P"' ■ (94) 
^ « ' ^ . ' 

For example, for x = — Ar this is just the Hamiltonian of the supersymmetrized Hydrogen 
atom which has been introduced and solved in It is evident from the representations 
and (|M)) that Q decreases and Q"^ increases the eigenvalue of the number operator 

iV = ^^"t^n (95) 

by one and H commutes with A^. The eigenvalues of N are 0, 1, ... , 2N. 

As before, we interpret the 2N coordinates x"" and annihilation operators ip" as values 
of two scalar and one Dirac field on a one-dimensional lattice with lattice sites. More 
precisely, we make the following identifications for n = 1,2, . . . , N , 



Hn) = f^^^l') , i^\n) = (^^^Z') ■ (96) 



The free supercharge H93() takes the form 



n=l ^ ' n=l 



The remaining task is to find a superpotential x giving rise to interacting lattice Wess- 
Zumino models. Since x should be real we use a representation for the two-dimensional 
7-matrices such that 17* and 7*^ are real, 

7° = o's, 7^ = 7* = -<72, (98) 

in order to obtain a real Dirac-Hamiltonian, 

As explained above, d need not be anti-Hermitian in which case we take 

hT = = -i7*9A + ^7° - ij^ds with (/i^)2 = (-A + m^)l2, (99) 

such that is real and Hermitian. Note that the term containing is not a momentum 
dependent mass term as in ()39() . We have been lead to a different type of Wilson term 
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as compared to the N = \ model since we have chosen a different representation for 
the 7-matrices. The earlier Majorana representation (|17() is not useful in the present 
context, since it would lead to a complex x iii (ESI)- 

The term in (|94|) must contain the free Dirac-Hamiltonian and this condition implies 

^2 777, 

- 75:rT^lxT7T = «'/5=1'2, n,n' = l,2,...,iV. (100) 

d(pa{n)d(j)p{n') 

Hence we expect that the real function 

x'^ = -\{ct>^K4>), (101) 

is the superpotential for a M = 2 supersymmetric model. For these models we use the 
following inner products 

{4>,4>) =^4>a,n4>a,n and (V', V') = ^ ("-) V'a ("-) , (102) 
a,n a,n 

for scalar doublets and Dirac spinors, respectively. The corresponding supercharge 

n \ / 

and its adjoint give rise to the following super Hamiltonian, 

= -^(^,vr) + i((/>,(-A+m2)0), H^ = {i;,h^ij). (104) 

Note that the superpotential is a harmonic function, Ax™ = 0, and thus there is no 
constant contribution to H^. The charges and Hamiltonian act on function(al)s in the 
Hilbert space of the M = 2 lattice models 

7^ = f)®...0(), where t) = LiiR^) (S) C"^ (105) 

^ V ' 

Af— times 

is the Hilbert space for the degrees of freedom on one lattice site. 
Now we turn to interacting models by replacing the mass term in 

n 

given by the quadratic harmonic function /, by an arbitrary harmonic function /(0) of 
the two variables (pi and (j)2, 

X = -U^,hU) + J2fi^(ri)), where A/ = 0. (106) 

The supercharge and its adjoint are calculated as 

Q = e-'^Qoe^ and = e^Qle-^ (107) 
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with Qq and qJ from 1)9 7|) . After some algebra one finds for the bosonic part of H 
^{(5,(5^} the formula 



^-^2) (108) 



where the harmonic functions / and g are the real and imaginary parts of an analytic 
function, such that 

df dg df dg 

0(pi a(p2 0(j)2 0(j)i 

For the fermionic part of the Hamiltonian one obtains 

/7p = (V,/i°V')-(V',7°r(0)^), r(</.) =/,n(0)- i7*/,i2 (</>). (110) 

The last term contains the Yukawa coupling between scalar and Dirac fields. Note that 
the last two terms in (|l()8j) can be rewritten as 

'^,d.A + i'^,^.d,A. (Ill) 



In the continuum limit the first term on the right becomes a surface term commuting 
with the supercharges and the second term, which is a lattice artifact, must vanish. Thus 
it is natural to set 

H^ + Hp = Pq + Z (112) 

and interpret the first term 

Po = i(7r, vr) - i(<A, A<P) + (V^, h^) + ^ (||' ^) " (^'^ 7°r(,/.)V) (113) 

as energy and the second term as 'would be' central charge Z in 1)111(1 . This agrees with 
our interpretation for solitonic configurations saturating the BPS-bound. To see that 
more clearly we consider the energy of a purely bosonic static solution, 

E = P„ = -i(*.A*) + ^(|.|)- (114) 

From the very construction it is evident, that there is a BPS-bound. One just adds the 
non-negative operator in H108|l to the non-negative operator one gets when changing 
the signs of / and g and finds 

E>\Z\. (115) 
For example, a cubic superpotential f + ig = A0^/3 leads to a (;/)'^-models with 

Po = ^(vr,^) - i(0, A(/<) + {^,hU) + i4>Af - (V',7°r(0)V') . (116) 
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It contains a scalar and pseudoscalar Yukawa interaction with 

r((/)) = 2A((/>i +i7*02). 
The would-be central charge is cubic in the scalar fields and reads 

Z = 2\ ((/>i(/.2, - A - <l)l d<t>2) ■ (117) 

Before turning to the discussion of the ground state we note, that the conserved number 
operator 

N = ^ij\n)^{n) (118) 

n 

leads to a decomposition of the Hilbert space (|1U5|1 in orthogonal subspaces labelled by 
the fermion number, 

n = no(Bni®---®n2N~i®n2N, n\ =pi. (ii9) 



The nilpotent supercharge Q decreases N by one, increases it by one and the super- 
Hamiltonian commutes with A^, 



[N,Q] = -Q, [iV,Qt]=QT and [N,H]=0. (120) 

We call the subspace Tip p-particle sector. The states in the zero-particle sector are 
annihilated by Q and those in the 2A''-particle sector by . 



5.2 Ground State of the Free Model 

The Hermitian lattice Dirac-Hamiltonian in (|99j) is real and can be diagonalized by 
an orthogonal matrix S, 

K^ = S-^DS, D = dmg{di,d2,...,d2N)- (121) 

We rotate the field-variables with S, 

^ = S<j), rj = Sip and r]^ = Sip^ . 

The new fields still obey the standard anticommutation relations, e.g. 

{r]l{n),r]f^{m)} = SafsSnm, (122) 

and the transformed supercharges read 

Q = i??(^-^e) and Qt = i^t|^^ + 2^^^ (123) 
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and show, that the new degrees of freedom decouple. Hence the ground state must have 
the product form 



*o = exp(-i^|d„|e^)|0) (124) 
and the supercharges act on this state as fohows, 

Q^o = 2i daVa^a^o , Q^^'o = -2i J] ^r/^a^o. (125) 

a:da>0 a:da<0 

This way we arrive at the following conditions for this state to be invariant, 

da>0^ r]a\n) = and da < ^ r]l\n) = 0. (126) 
This leads to the unique normalizable ground state (|124jl with 

\^) = n (127) 

da<0 

which is annihilated by the supercharges and hence has vanishing energy. There are N 
positive and N negative eigenvalues of such that the invariant vacuum state lies in 
the middle sector Ti^ in the decomposition ()119() of the Hilbert space. All fermionic 
states with negative energies are filled. This is just the Dirac-sea filling prescription. 
Note that our result is the lattice version of the continuum result for the ground state. 



^'o = exp (^-^ j (t>aV-^ + m'^'Pa^ 



5.3 Ground States for Strong Coupling 

In the strong-coupling limit we may neglect the spatial derivatives such that the super- 
charges and the Hamiltonian becomes the sum of N commuting operators, each defined 
on one lattice site [2j|. The operators on a given site take the form 

Q = iV,(V + V/), i/ = -iA + i(V/,V/)-VV'V. (128) 
Now we explicitly construct the ground state for the harmonic superpotential 

/(0) = -Re(?^, <^ = 01 + i02 = re'^ (129) 
P 

which gives rise to a super symmetric anharmonic oscillator on the Hilbert space f) = 
L2(1R.^) X C^. The bosonic part of H reads 

Hs = -^A + V with V = iAV2p-2, (130) 
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and its fermionic part 



n \( U /t /^-Re(/>P-2 Imr-^\ 

H. = X{p-l)i^'[^j^^-2 ^,^^2)i^ (131) 

It is useful to note that H commutes with the operator 

J = L + S, S = -sip^a2%^, s = i(p-2) (132) 

and that the ground state must reside in the two-dimensional sector with particle number 
since the restriction of Hp to the zero- and two-particle sectors vanish 
and i^B > 0. The one-particle sector is spanned by the following two eigenstates of S, 

I T) = (V'l - iV'l) |0) and I i) = (V't + i^t^ |o) (133) 

with eigenvalues 1 and —1, respectively. Here |0) denotes the Fock- vacuum which is 
annihilated by the annihilation operators ipa- Diagonalising J in this sector leads to the 
ansatz 

Vo,(0) = Rj+{r)e'^^-'^' I T) + i?,_(r)e^(-'-+^)^ | j), (134) 

where the J-eigenvalue j is integer for even p and half-integer for odd p. Inserting into 
Qijjj = Q^ipj = yields the following coupled system of first order differential equations 
for the radial functions 

i?;.±(r) - '-^Rj±{r) + Xr'-'R,^{r) = 0. 

The square integrable solutions are Bessel functions 

Rj±(.r) =crP-^K,^^(^^rP') with j e {-s,-s + 1, . . . ,s - l,s}. (135) 

The number of supersymmetric ground states of the models with (fP''P~'^ self-interaction 
is just p — 1. The {p — 1)^ normalizable invariant eigenstates are 

N 

^0,h,-dN = (8)'^0i„(<^n) G f)l X • • • X 1)^. (136) 
n=l 

For example, for the M = 2 model with cf)'^ interaction there exist 2-^ normalizable zero 
modes in the strong-coupling limit. This number diverges in the thermodynamic limit. 
On the other hand, there is exactly one normalisable zero mode when one switches off 
the self-interaction. This discrepancy between the number of supersymmetric ground 
states in the weak and strong-coupling regimes becomes even more puzzling when one 
takes into account certain rigorous theorems on the stability of such states under analytic 
perturbations discussed in the following section. The zero modes in (|136|1 with radial 
functions (|135|) have been constructed previously in [221331 and [HK] . 
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6 From Strong to Weak Couplings 



6.1 Perturbation Theory and Zero Modes 

Let us recall a well known result for perturbation theory of zero modes in super symmetric 
quantum mechanics [33] • We consider the N = 1 case and denote the single Hermitian 
supercharge by Qq, 

Ql = Ho, {r,go} = o, rt = r, = i. (137) 

In addition, we define the projection operators 

V± = ^{1±T) (138) 

which project on the ±1 eigenspaces of F. These eigenspaces are denoted by TCb/f- Iii 
the following we assume that there are no zero modes in Tip and at least one zero mode 
^'o in the bosonic sector TCb- We perturb the operator Qo by an operator eQi with real 
parameter e, Q(e) = Qo + eQi, where {Qi,T} = 0. We want to solve the eigenvalue 
equation 

Q{e)^{e) = A(e)^(6), 
with A(0) = and ^^(O) = ^I'o- We consider the following formal power series in e, 

oo oo 
k=l k=l 

Proposition: Under the assumptions above one has A(e) = and r^'(e) = ^'(e) in the 
sense of formal power series. 

Proof by induction: To order e'^ the proposition holds. Let us assume it holds up to 
order e^~^. To order we obtain the equation 

Qo^i + Qi*i-i = A,^o. (139) 

Taking the scalar product with ^'o yields 

Xj = (^'o,Qi*j-i). 

Since T squares to 1 and anticommutes with the perturbation we find 

A, = (r2^'o,Qi^',-i) = -(r^o,Qir^',-i) = -(^cQi^i-i) = -A,-, 

which proves that Xj = 0. Furthermore, with 
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we conclude 



QoV-^j = 0, 

where we used the projection operator V- introduced in ()138() . As V^^j is a zero mode 
of Qo it fohows by assumption that it resides in TCb- But as V- projects onto Tip, 
we conclude V-'^j = or E TC^- This then proves our statement. Note that the 
statement has been proved in the sense of formal power series only. In case A(e) is not 
analytic at e = the result above maybe misleading. 

6.2 The X = 1 Case 

In what follows, we compare the strong-coupling results with the usual perturbation 
theory around minima of the potential. 

In the case deg(iy) = p even, supersymmetry is never broken, neither in the strong- 
coupling limit nor in perturbation theory. For even p there is at least one minimum of 
the potential V = ^{W')"^ with V = 0. The quadratic approximation of the potential at 
the critical points yields for each minimum one normalizable zero mode similar to the 
ground state of the free model. In contrast to the strong-coupling limit there may be 
more than one perturbative zero mode, but they always come in an odd number. The 
difference of bosonic and fermionic zero modes is ±1 as in the strong-coupling limit. 

In the case deg{W) = p odd, the difference between the strong-coupling limit and per- 
turbation theory is more severe. Supersymmetry is broken in the strong-coupling limit 
but it may be unbroken in perturbation theory. Let us consider an explicit example, 

W{<P) = + goci). (140) 

Perturbation theory for go < predicts one bosonic and one fermionic zero mode (un- 
broken supersymmetry), and broken supersymmetry for qq > 0. The strong-coupling 
limit states that supersymmetry is broken for all qq. 

In Appendix IB. II we provide the rigorous proof that \Ai with Ai given in H69|) is an 
analytic perturbation of Aq in (|69|) . This implies that all eigenvalues are analytic func- 
tions of the parameter A. Assume now that in a finite range of the parameter A there 
is a ground state with energy zero. As an analytic function which vanishes in some 
finite range is identically zero, the number of zero modes changes at most at isolated 
points of the parameter space of A. Furthermore, in the strong-coupling limit, we have 
either bosonic or fermionic zero modes. In subsection 16 . 1 1 we have proved that under this 
assumption a zero mode always remains a zero mode. We conclude that, generically, 
the number of zero modes is given by the number of zero modes in the strong-coupling 
limit. Generically, since for certain discrete values of A the number of zero modes could 
be enhanced. Moreover, as the index also depends analytically on the parameter A, we 
are able to calculate this index in the strong-coupling limit. 
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In the continuum and infinite- volume limit these arguments may break down, as the 
estimates necessary for proving analyticity (see Appendix IB.lj) may not be valid anymore. 
In the unbroken case we can definitely conclude that the theory is still unbroken in the 
continuum and infinite- volume limit. Suppose we know that for any finite lattice there 
is at least one ground state with zero energy. As the limit of zero is again zero this 
mode survives in the limit. In the case of broken supersymmetry a non-zero energy 
eigenstate may become a zero mode in the continuum and infinite-volume limit, and 
supersymmetry may get restored in this limit although it is broken for all finite lattices. 

Indeed, for negative go in our example above, the scalar field has a non-vanishing vacuum 
expectation value and therefore the fermionic field tp acquires a non-zero mass. As there 
is no massless Goldstone fermion, supersymmetry has to be unbroken in this case 

Let us summarize. On a finite lattice, the strong- coupling limit gives the correct number 
of zero modes of the full problem. There is only one zero mode in the case where 
deg(VF) = p is even, and otherwise there is no zero mode. Variations of the parameters 
in the superpotential of power less than p do neither change the number of zero modes 
nor the index. For example, in the model with superpotential given in H14U|) . it is 
impossible to have two phases of broken and unbroken supersymmetry (depending on 
the value of the parameter go) on a finite lattice. The numerical simulations in |^ may 
be interpreted as hinting towards such a phase transition in the continuum theory. 

6.3 The M =2 Case 

Similar to the case = 1, we prove in Appendix IB. 21 that the index in the strong- 
coupling limit is the same as for the full problem. This implies that we have at least 
{p — 1)^ zero modes for the theory on finite lattices. For the continuum theory in a 
finite volume, it was shown using methods of constructive field theory that the N = 2 
Wess-Zumino model is ultraviolet finite and that the index is given by p — 1 |24j . This 
seems to be in contradiction with our result, as the {jp — 1)^ zero modes exist for all 
finite lattices and, by the same arguments as for the N = \ model, remain zero modes 
in the continuum limit. 

We suggest the following solution for this problem. Remember that our lattice Hamilto- 
nian H contains not only the discretized version of the continuum Hamiltonian Pq but 
also the central charge Z, i.e. 

H = Po + Z. (141) 

Furthermore, both Pq and Z contain the lattice derivative which couples fields at different 
lattice sites. If we choose in the strong-coupling limit a zero mode that varies from lattice 
point to lattice point, both Pq and Z may become very large but will, nevertheless, add 
up to zero. In the continuum limit the energy Pq may be infinite in which case this 
rapidly varying zero mode is only a lattice artifact. On the other hand, if we choose 
the same zero mode for each lattice site, then Pq as well as Z should be zero in the 



28 



continuum limit. Thus, there are exactly p — 1 such modes. We are planning to test 
this proposal in a perturbative calculation of Pq and Z. The results will be presented 
elsewhere. 

7 Conclusions 

In this paper we have related Dirac operators defining supersymmetric quantum me- 
chanical systems in high-dimensional spaces [32, to Wess-Zumino models on a spatial 
lattice in 1 + 1 dimensions. After a very particular dimensional reduction the square 
of ilp can be identified with the super-Hamiltonians of latticized Wess-Zumino models. 
This way we discovered a natural connection between discretized supersymmetric field 
theories and supersymmetric quantum mechanics. 

We have recalled the continuum formulation of Wess-Zumino models and discussed their 
lattice versions. For the case of simple {M = 1) and extended (AA = 2) supersymmetry, 
we have derived the corresponding Dirac operators. Furthermore, all ground states for 
the free massive models and the interacting theories in the limit of strong coupling have 
been constructed. 

Different realizations of lattice derivatives have been discussed and their properties - in 
particular from the point of view of supersymmetric quantum mechanics - have been 
analyzed. Our results on the number of zero modes do not depend on the particular 
lattice derivative, as long as some mild assumptions are fulfilled. 

Employing powerful theorems from functional analysis we were able to relate the strong 
and weak coupling regions. For = 1 it turns out that generically the number of zero 
modes is determined by the strong-coupling limit. There is a single (no) zero mode, if 
the degree of the superpotential is even (odd). For M = 2 we find at least {p — 1)^ zero 
modes, where p is the degree of the superpotential and N the number of lattice sites. 
This number is far off the correct continuum result, which predicts p — 1 zero modes, a 
serious problem which has been observed earlier in j21j . 

We have explained this paradox as follows: the lattice Hamiltonian H does not only con- 
tain the continuum Hamiltonian Pq but also additional terms which (for antisymmetric 
lattice derivatives) are to be interpreted as a lattice version of the central charge Z. On 
the lattice, Pq and Z cancel pairwise for the huge number of zero modes under discus- 
sion, even though neither Pq nor Z is zero in the continuum limit, except for exactly 
p — 1 of the modes. 

Our Dirac operators clearly deserve further studies. For instance, the application of 
(standard) index theorems to the case at hand should reveal new information about the 
field theories. We also plan to extend our results to Dirac operators on curved manifolds, 
which can be reinterpreted as nonlinear cr-models from the field theory point of view. 
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A The SLAC Operator 

In this appendix we introduce and discuss the nonlocal SLAC derivative. It can be used 
to define chiral fermions without fermion-doubling. 

First we consider real valued scalar fields on the spatial lattice. They maybe interpreted 
as wave functions of a quantum mechanical system with Hilbert space IR^, equipped 
with the scalar product 

N 
n=l 

Although the fields are real it is useful to embed them in the space of complex valued 
lattice fields. For a normalised function we interpret |(/>(n)p as probability for finding the 
'particle' at the lattice site n. Expectation values of functions of the 'position' operator 
n are 

{f{n))^ = Y,Hn)fin)4>{n). (142) 

We want to derive a similar formula for expectation values of functions of the momentum 
operator. For that aim we Fourier transform the wave function as follows 

N' 

4>{p,) = ^ ^ e-'f'=-</.(n), where AT' = i(Ar-l), Pk = ^k. (143) 

'^^ n=-N' 

The inverse Fourier transformation reads 

1 

•^^ = 7^ E ^""-""^iPk)^ n = -N',...,N'. (144) 

V ^ k=-N' 

We have chosen the symmetric summation to end up with a real difference operator. For 
periodic fields n must be integer and this is the case for space lattices with an odd number 
of sites. For a normalized the Fourier transform cj) is normalized as well and we may 
interpret |0(pjfc)P as probability for finding the 'momentum' p^. With this interpretation 
we obtain 

N' 

{f{p))<t>= Yl fiPk)\4>iPkt = Y^nfip)nn'4>n' (145) 
k=-N' nn' 
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with matrix elements 

N' 

N 



1 ^ 



k=-N' 

With the help of the generating function 

N' . , . 

Z{x)= V ^iP.^ = ^}p^ (146) 

we can calculate all matrix elements of f{p)- Now we are ready to define the real, 
nonlocal and antisymmetric lattice operator 9slac = ip- The matrix elements are 

/(4..c)..„. = i,/ (^) Z(.)l^^_^, . (147) 

As expected 9slac is a Toeplitz matrix with elements 

(^SLACW = (-)"-"' . . (1^^) 

sm (^7r(n — n')/jyj) 
and the elements on the diagonal vanish, (5sLAc)nn = 0. 



B Analyticity of Perturbation 

In the following we consider operators on the Hilbert space 

7^ = L2(R'^,d'^x)®C^ (149) 

for Dsn with norm 

D 

ll/f = Ell/^llL> /=(/i,...,/d)gW. (150) 

1=1 

Here, || • ||l2 denotes the familiar L2-norm. 



B.l The X =1 Case 

For the Wess-Zumino model on the lattice with Af = 1 supersymmetry we take D = 2^ , 
d = N {N = number of lattice points) and consider the (unperturbed) operator H69|) 

N 

^o = Yl {-^Mn)dn + Mn)W'{xn)) ■ (151) 

n=l 
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We recall that is a polynomial of degree deg(VK) = p > 1 and ipain) are Hermitian 
D X D— matrices obeying the Clifford algebra 

{Mn),4^^{n')} = 26a/3S{n,n'), a,/3 = l,2, n,n' = 1,2, N. (152) 

The operator Aq with domain of definition 

V{Ao) = C^iB.^) 

is essentially self-adjoint, where we write C^(IR^) for the C°°-functions with compact 
support in IR.^ . A simple calculation using (|152() shows 

(^O)' = Yl i-^r.dn + W'{Xn)W'{Xn) " i^l (n) ^2 (n)Ty"(x„)) . 



Closure of the Operator 

To determine the closure of the operator Aq we have to find the closure of its domain 
^{Aq) with respect to the norm 

\\f\\M,a = a\\ff + \\Aoff, a>0. (153) 

Note that these norms are equivalent for all a > 0. Using the abbreviation 

Pp = l + \x\P-\ (154) 

we can prove the following 

Lemma: There exist constants a, 6i, 62 > such that 

Wf'f + hWPpff < ll/llio,. < Wf'f + b2\\Ppff (155) 
holds for all f £V{Ao). 

In the Lemma we used the short hand notation = Ylm 

Proof: First, we show that only the degree deg(Ty) = p is important for terms like 
En W{xn)ff with / E V{Ao). Indeed, we find 



J2\\W{xn)ff<Na,\\ppff, a. 



W'iXn) 



Pp 



(156) 



where || • ||oo denotes the supremum norm. The factor N arises from the sum over n, as 
oi does not depend on n. Similar we obtain 



\Ppft 



Pp 



l + Y,W'^Xn)f 



v^+j:„w'Hxn. 

< a2(\\ff + Y,\\W'{x^,)ff), 02 



Pp 



■ (157) 
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Now, it is easy to prove the second inequality in (|155() . 

ri5fil ^ 

a||/f + \\Aoff < WfT + + Na,\\p,ff + ^ ||/|| \\W"{x^)f\\ 



< 



P'l|2 



+ {a + Nai+Na3)\\ppf\\' 



(158) 



62 



with 03 



values are ±1. In t 



. We used that the matrix- norm of ipa{n-) is one, since its eigen- 
le last inequality we made use of ||/|| < which holds for all 



The other inequality in H155|) is more difficult to prove. With (|157|) we get 

+ \\Aoff > Wf'f + -WPpff + (a - - 11/11 \\W"{xn)f\\. (159) 
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In order to obtain a positive constant bi in our lemma we must be rather careful with 
our estimates for the last term in (|159j) . We introduce a ball of radius R and split / into 
two parts, f = f< + f>, where /< has its support inside the ball /> outside the ball. We 
obtain 



11/11 ||Ty"(x„)/|| = J2 (ll/<ll ll^^"(^n)/<|| + ||/>|| \\W"{x„)f: 



>\\ > 



(160) 



where the terms containing both /< and /> vanishes. Let us now consider the two terms 
separately. First, we obtain 



Yl \\f>\\ ll^"(^n)/>|| < Nai{R)\\ppff, a4(i?) 



W"{Xn) 



Pp 



(161) 



00, > 



where we have introduced the supremum norm || • ||oo,> = sup|^|>^{| • |}. For large R we 
have 04 (i?) ~ 1/R such that 04 gets arbitrarily small for big balls. Second, we obtain 



El 



|^"(x„)/<|| < Na^iRmf, a^iR) = \\W"{xn)\\oo,<, 



(162) 



with 



3o,< = sup|2,|<^{| • |}. For R ^ 00 we have a5(i?) — > 00. Altogether, we find 

+ \\Aoff > Wf'f + (a - 1 - Na,iR)) \\ff + (1^ - Na,{R)) \\ppff. (163) 

^ V ' 

61 

In a first step we choose R large enough such that bi is positive. In a second step we 
choose a large such that the constant in front of ||/|p is positive as well. This finishes 
the proof of our lemma. 

Since all norms 



(164) 
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are equivalent for 6 > 0, the Lemma implies that these norms are equivalent to the 
norms „ in H153|l . Therefore, the closure of P(^o) with respect to the norm (|153() 

coincides with the closure with respect to || • ||fe. This closure is given by 

P(Ao) = {/ G WiiR"") : IIpp/II < oo} ^ W^iR'',pl) C^. (165) 

Here W^CR^) is the Sobolev space with first weak-derivative in L2. 
Perturbation 

Let us perturb the operator ^0 by the operator Ai in (jSHI), 

N 

Ai = - Yl ^-(5L„V'2(n). (166) 

m,n=l 

The operator Ai is self-adjoint with 2?(^i) = L2(1R^,/5)®C^ D P(i^o), with p -weighted 
Lebesgue measure, where p{x) = (1 + |x|)^. From the following Lemma we will derive 
useful information about the nature of the perturbation. 

Lemma: For all A G R, and arbitrarily small e > there exists a > such that 

||A^i/|| <ePo/||+a||/||, WfGViAo). (167) 

Proof: We prove the inequality for all / G 'D{Aq). Then it holds for all elements in the 
closure as well. As before we split f = f< + f>- First, we note 

||AAi/<|| < \X\N\{R) 11/11, a{R) = \\xn\\oo,< ■ max{|a„„| : m,n = 1, . . . ,N}. (168) 

For R ^ 00 the constant a{R) tends to infinity. Next, we have 

||A^i/>|| < \X\N%{R){c\\f\\ + \\Aof\\), 



b{R) 



Pp{x) 



max{|Omn| : n = 1, . . . , A^} (169) 



00, > 



for some positive constant c. For big R the constant b{R) tends to zero. We choose the 
ball big enoug h such that \X\N^b{R) = e and set Q = ec + \X\N^a{R). Note that the 
latter constant may become huge. 

Self-adjointness 

We did prove that Aq is a self-adjoint operator. Clearly XAi is symmetric on 'D{Ao). 
Furthermore, ()167() shows that A^i is Ag-bounded with relative bound less than one. 
The famous Kato-Rellich Theorem, see Theorem X.12 in |36j . states that under these 
conditions the operator 

Qi{X)=Ao + XAi (170) 

is self-adjoint with domain 'D{Aq). We conclude that Qi{X) is a family of self-adjoint 
operators with common domain of definition 'D{Aq). 
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Analyticity of Eigenvalues 

In the following we prove that Qi(A) is an analytic family in the sense of Kato for all 
real A. We have seen that Qi{X) is self-adjoint for real A. For a self-adjoint and analytic 
family it is known that the eigenvalues depend analytically on the parameter A, see for 
example Theorem X11.13 in |36| . 

For an arbitrary real Aq the perturbation Aq^i is ^o-bounded with arbitrary small 
relative bound ()167|) . Then, it is easy to see that Ai is (5i(Ao) -bounded. It follows 
that for small e the operators Qi{Xq + e) form an analytic family of type (A) [HE] and 
therefore also an analytic family in the sense of Kato. But as Aq G K. is arbitrary, we 
haven proved that Qi{X) is analytic for all real A. 

Actually, the cited Theorem X11.13 [^2] above is only valid for isolated eigenvalues with 
finite degeneracy or equivalently for eigenvalue in the discrete spectrum. In the following 
we prove that the spectrum of Qi{X) is discrete by proving this statement for its square, 
H{X) = (5i(A)^. H{X) is self-adjoint with domain of definition given by 

P(if(A)) ^ {/GP(Ao) :Qi(A)/GP(^)} 

= VF|(R^,p')®<D^, p'ix) = {l + \x\'P-')' (171) 

and it is semibounded 

H{X) > 0. (172) 

Such operators possess entirely discrete spectra if and only if its resolvent is a compact 
operator, see Theorem XIII. 64 in (36) . In the following we prove that i?(A) has compact 
resolvent for all A G R. 

We must prove that the image of a bounded subset of the Hilbert space, say 

{/GW: 11/11 <1}, (173) 

is mapped to a precompact set under the map {H — z)~^ for some z in the resolvent of 
H. The image is given by 

{g€V{H):\\{H-z)g\\<l}. (174) 

As earlier we split g into and g< and obtain for large enough radii R the inequality 
Ibll < lb<ll + €• Foi' a compact ball B = {x ^ : \x\ < R} we have Sobolev's 
embedding theorem and there is an e-net gj G (i^, p'), j = I, . . . , with \\g<c—9j II < e 
for one j G {1, . . . , A^}. We extend the gj by zero to the region outside the ball and 
obtain 

\\g-gj\\ <2e (175) 

for any g in the image of the unit ball under (H — z)^^ and a specific j G {1, . . . , N^}. 
We conclude that there is a 2e-net of the image and therefore the image is precompact. 
This completes our proof. 



35 



Stability of the Index 

We have shown that the eigenvalues are analytic functions of the parameter A on the 
whole real axis. It follows at onces that the index - the difference of bosonic zero modes 
and fermionic zero modes - is also an analytic function and, as the index only takes on 
integer values, is constant. 

An alternative and elegant proof of this statement can be given with the help of the 
theorem that a relatively compact perturbation does not change the index Indeed, 
inequality (|167|) implies that our perturbation is relatively compact^. 



B.2 The X = 2 Case 



After the detailed investigation of the = 1 case we shorten our discussion for J\f = 2. 
In what follows we consider the real part of the complex supercharge in H1U7() 

^0 = E f-iV'l(«)9x„ - i^Pl{n)dy„ + ^ V'2"(n)W^,a(Xn,yn)') (176) 
n=l \ a=l ) 

in the strong-coupling limit. For M = 2 supersymmetry the il)'^ are D-dimensional 
Hermitian matrices obeying the Clifford algebra, with D = 2^^ . The function W{x,y) 
is harmonic and and thus is the real part of an analytic function F{x + \y). As in chapter 
|2lwe use the notation VF,i (x,y) = dxW{x^y) and W,2 ix,y) = dyW{x,y). As domain of 
definition we take 

V{Ao) = C7~(R2iV) C^, D = 2^^, (177) 
such that Bq is essentially self-adjoint. We introduce the potential 

Kix,y) = Y,W,lix,y). (178) 

a 

For large radii r only the leading power of W is relevant. Therefore, we may consider 
the particular case 

W(x,y) = -RezP (179) 
p 

for which we find K{x, y) = KrP~^ ^ oo in all directions for r oo. 
The perturbation contains the lattice derivative, 

2N 



Bi= ^ (x.m{d)mn^2{'ri) + ym{0^)mni^l{n)] . 



m,n=l 

Replacing in the estimates for the case with M =1 supersymmetry the potential W'{xn) 
by K{xn,yn) leads to analogous results in the J\f = 2 case. Again all eigenvalues are 
analytic functions of the parameter A and in particular the index does not depend on 
this parameter. 

^We thank H. Triebel for the proof of this statement. 
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